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We study neutron matter at and near the unitary limit using a low-momentum ring diagram ap- 
proach. By slightly tuning the meson-exchange CD-Bonn potential, neutron-neutron potentials with 
various 1 5o scattering lengths such as a s = — 12070/m and +21 fm are constructed. Such potentials 
are renormalized with rigorous procedures to give the corresponding a s -equivalent low-momentum 
potentials Vi ow -k , with which the low-momentum particle-particle hole-hole ring diagrams are 
summed up to all orders, giving the ground state energy Eo of neutron matter for various scattering 
lengths. At the limit of a 3 — > ±oo, our calculated ratio of Eo to that of the non-interacting case 
is found remarkably close to a constant of 0.44 over a wide range of Fermi-momenta. This result 
reveals an universality that is well consistent with the recent experimental and Monte-Carlo compu- 
tational study on low-density cold Fermi gas at the unitary limit. The overall behavior of this ratio 
obtained with various scattering lengths is presented and discussed. Ring-diagram results obtained 
with Viom-k and those with G-matrix interactions are compared. 
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I. INTRODUCTION 

Back in 1999, BertschQ formulated a many-body 
problem, asking: what are the ground state properties 
of a two-species fermion system that has a zero-range in- 
teraction and an infinite scattering length? Such problem 
was originally set up as a parameter-free model for a fic- 
titious neutron matter. Recently, as the experiments on 
trapped cold alkali gas undergo huge breakthroughs, de- 
generate Fermi gas with a tunable scattering length (in- 
cluding ±oo ) becomes accessible in laboratories [2|. Since 
then cold Fermi systems have aroused growing attention. 

The term 'unitary limit' has been used by many au- 
thors to refer to the special scenario in a low-density two- 
species many-body system where the scattering length 
between particles approaches infinity. More specifically, 
at the unitary limit, the scattering length a s , the Fermi 
momentum kp, and the range of the interaction r m + sat- 



isfy a J >> k F » r: 



Under such condition, atoms 



are 'strongly interacting', and a full theoretical descrip- 
tion of their properties is a challenging task in many-body 
theory. Universal behavior is expected to show up in 
various aspects, including ground state properties as dis- 
cussed below, collective excitations H, EIlE la ^ IE M , and 
thermodynamic properties [l(| HH, IeOiS |24 1 . Such uni- 
versality can be naively understood as the 'dropping' of 
the scattering length o s out of the problem, leaving kp as 
the only relevant length scale. In particular, the ground 
state energy Eo, is expected to be proportional to that of 
the non-interacting gas Eq 
, or equivalently 

Eo 
A 



free 



[3], that is E /E f 



free 



3 k F 
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(H = m = 1), A being the number of particles. The uni- 
versal constant £ is of great interest and many attempts 
have been made to derive it analytically or determine it 
exp er iment ally. 

Theoretical calculations suggest that £ is between 0.3 
to 0.7. For example, an early work based on dif- 
ferent Pade approximations gives £ = 0.326. 0.568[14j. 
Diagrammatic approach gives 0.326 with Galitskii 
resummationfTH . 0.7 with ladder approximation [l5j], 
and 0.455 with a diagrammatic BCS-BEC crossover 
theory [l7j. Other theoretical approaches have also been 
used, including e expansion, which gives £=0.475 in 18] 
and [20|| , and variational formalism, which gives 0.360 in 
[l||. The four most recent experimental measurements 
are listed in Table |TJ Though the experimental results 
are consistent with each other, the experimentally deter- 
mined value of £ still falls between relatively large error 
bars(~10%). By far the best estimate on £ is considered 
to be that from Quantum Monte-Carlo methods, giving 
£ = 0.44(1) [Ij and 0.42(1) [l^. 



£ 



Authors 



Ref. 



0.36(15) 
0.51(4) 
0.46(5) 

0.46 



0.05 
0.12 



Bourdel el.al 
Kinast et.al. 
Partridge et.al. 
Stewart et.al. 



[23] 
[24] 
f25] 
[26] 
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TABLE I: Comparison of recent experimental values on £. 

Cold and dilute neutron matter is a special class of 
cold Fermi system with great importance in astrophysics. 
Its properties at resonance has attracted much inter- 
est recently [13, H3- ^ n tn i s work we report results 
from low-momentum ring diagram calculations on the 
ground-state energy of neutron matter at and near the 
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unitary limit. As is well-known, the 1 So channel of 
neutron matter has a fairly large scattering length a s 
(— 18.97/m), nonetheless, it is still finite. Here, by ad- 
justing the interaction parameters of the CD-Bonn po- 
tential [29( | . we construct 'tuned' neutron interactions 
with different a s 's such as — 9.83/m, — 12070/m and 
+21 fm (which possesses a bound state). For a wide 
range of neutron density, the case of a s = — 12070/771 
can be considered the same as the unitary limit, namely 
a s — > — oo. We shall compute the ground state en- 
ergy of neutron matter, with inter- neutron potentials 
being these 'tuned' CD-Bonn's, by two steps: renor- 
malization followed by ring summation. We first renor- 
malize neutron interactions with a T-matrix equivalence 
renormalization method [3(J HH, [H, [3(1], where 

the high-momentum components beyond a decimation 
scale A are integrated out. This gives the corresponding 
low-momentum interactions Vi ow - k s with the scattering 
lengths being preserved. Then, we calculate the ground 
state energy by summing the particle-particle-hole-hole 
(pphh) ring diagrams]!?} to all orders. In such ring sum- 
mation, we employ a model space approach, namely, the 
summation is carried out within a model space charac- 
terized by {k < A}. 

We shall closely examine how our results differ from 
similar calculations with a different renormalizcd inter- 
action - the Brueckner G-matrix on which the Brueckner 
Hartree-Fock(BHF) method is based. The BHF method 
has been widely used for treating the strongly interacting 
nuclear many body problems (38.|39|. However, BHF is a 
lowest-order reaction matrix (G-matrix) theory and may 
be improved in several aspects. To take care of the short 
range correlations, the ladder diagrams of two particles 
interacting with the bare interaction are summed to all 
orders in BHF. However, this method does not include 
diagrams representing hole-hole correlations such as dia- 
gram (iii) of Fig. 1 . Note that this diagram has repeated 
(pphh) interactions as well as self-energy insertions to 
both hole and particle lines. Another aspect of the tra- 
ditional BHF is that it employs a discontinuous single- 
particle (s.p.) spectrum which has a gap at the Fermi 
surface hp. To improve upon these drawbacks, Song et 
al. [371 ] have formulated a G-matrix ring-diagram method 
for nuclear matter, with which the pphh ring diagrams 
such as diagrams (i) to (iii) of Fig.l are summed to all 
orders. This ring-diagram method has been applied to 
nuclear matter and given satisfactory result [37j. The 
Viow~k ring diagram method used in this work is highly 
similar to [371 ] 's , except for one significant difference: the 
interaction used in the G-matrix ring diagram method is 
energy dependent. (The Brueckner G-matrix is energy 
dependent, as we shall later discuss.) This complicates 
the calculation a lot. Vi ow - k provides a cleaner and sim- 
pler implementation on such all-order ring summation. 

We shall first provide an outline of the ring-diagram 
approach in section II. The derivation details of the 
low-momentum interaction from the CD-Bonn potentials 
shall be followed in section III. Our major results from 



the Viow-k ring diagram method are in section IV. There 
we shall present our results for the ground-state energy 
and ratio Eq/Eq 7 "^ obtained with potentials of various 
scattering lengths. A fixed-point criterion for determin- 
ing the decimation scale A will be discussed. There one 
can also find a comparison of data on the ground state 
energy obtained with two different methods-the Viow-k 
and the G-matrix ring diagram methods. We shall sum- 
marize and discuss our work in the last section. 



II. LOW-MOMENTUM RING DIAGRAMS 

In this section wc describe how we calculate the ring 
diagrams for the ground state energy shift AEVj, which 
is defined as the difference (Eq — EQ Tee ) where Eq is the 
true ground-state energy and E^ ree is the corresponding 
quantity for the non-interacting system. In the present 
work, we consider the pphh ring diagrams as shown in 
Fig. 1. We shall calculate the all-order sum, denoted 
as AEq P , of such diagrams. Our calculation is carried 
out within a low-momentum model space {k < A} and 
each vertex of the diagrams is the renormalized effec- 
tive interaction corresponding to this model space. Two 
types of such interactions will be employed, one being 
the energy-independent Vi ow -k and the other being the 
energy-dependent G-matrix interaction. Let us consider 
first the former. In this case, AE PP can be written [37| 
as 



AE PP 



, , , due™ tr <A [F(u)V lc 
+ ^(F(cj)V low - k ) 2 + ^{F(Lj)V low - k ) 3 + 



where F is the free pphh propagator 



F ab (u) 



n a n b 



n a n b 



lu - (e a + e b ) +i0+ uj - (e a + e b ) - i0+ 



(3) 



with n a = 1, a < hp; — 0, k > kp and n a = (1 — n a ). 

We now introduce a strength parameter A and a A- 
dependent Green function G pp (uj, A) defined by 

G pp (lu, A) - F(w) + \F(lu)V 1ow ^G pp (lu, A). (4) 

The energy shift then takes the following simple form 
when expressed in terms of G pp , namely 



r / / .i,.,n+ 



AE PP = — ' "i 
2tti Jo 



o - — / dX I e tr <A [G pp (u,\)V low - k ] 



Using Lehmann's representation for G pp , one can show 
that 



AE% P = / dXE m 'Eijki<AY m (ij, X)Y^(kl, X)(ij\Vi ow - k \kl), 
Jo 

(6) 
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where the transition amplitudes Y are given by the fol- 
lowing RPA equation: 

yitfa + £ j)Sij,ef + A(l - rii - nj)(ij\Vi ow - k \ef)] 

ef 

xY m (ef,X) =u m Y m (ij,\); (i,j,e,f)<A. (7) 

The index m denotes states dominated by hole-hole com- 
ponents, namely, states that satisfy (Y m \-^\Y m ) = — 1 
and Q(i,j) = (1 — rij — rij). We have used the HF s.p. 
spectrum given by Vz OU) -fc, namely 

e k = h 2 k 2 /2m+ J2 (kh\V low - k \kh) (8) 

h<k F 

for both holes and particles with k < A. Thus the propa- 
gators of the diagrams as shown in Fig. 1 all include HF 
insertions to all orders. The above spectrum is continu- 
ous up to A. 

The above ring-diagram method is a renormalization 
group approach for a momentum model space defined by 
a momentum boundary A, and the space with momentum 
greater than A is integrated out. The resulting effective 
interaction for the model space is Vi ow ~k which is en- 
ergy independent. This renormalization procedure can, 
however, also lead to a model-space effective interaction 
which is en erg y dependent. The G-matrix ring-diagram 
method of [37J is of the latter approach. Formally, these 
two approaches should be the same. In the present work 
we shall carry out ring-diagram calculations using both 
approaches; it would be of interest to compare the results 
of these two different approaches. 

In the following, let us briefly describe the G-matrix 
ring diagram method [37| ■ Here each vertex of Fig. 1 is a 
model-space G-matrix interaction, to be denoted as G M . 
It is defined by 

Gfjkliu) = KjM + V vrs ~ _ p Jfl\ iQ + G rlkl M 

rs r s 

(9) 

where k 2 stands for the kinetic energy h 2 k 2 /2m and sim- 
ilarly for kg. The Pauli projection operator Q M is to 
assure the intermediate states being outside A and kp, 
namely it is defined by 

Q M (rs) = I, if max(k ri k s ) > A and min(k r , k s ) < kp 
= 0, otherwise. (10) 

In the above kp < A. In Rcf.[37] A is chosen to be 
~ 3/m -1 . Note that the above G M is energy dependent, 
namely it is dependent on the energy variable to. How- 
ever, ui is not a free parameter; it is to be determined 
in a self-consistent way. For example, the model-space 
s.p. spectrum is given by the following self-consistent 
equations: 

^ = ^3 + {a\U\a); (11) 



(a\U\a) = (a,h\G M (u = e a + e h )\a,h), a< A 

h<kp 

= 0, otherwise. (12) 

In the above U is the s.p. potential and e the model-space 
s.p. energy which is determined self-consistently with the 
energy variable of G M . Note that this s.p. spectrum 
does not have a gap at kp; it is a continuous one up to 
A. When choosing A=kp the above is the same as the 
self-consistent BHF s.p. spectrum. 

When calculating the ring diagrams using G M , its 
energy variable is also determined self-consistently. In 
terms of G M , the all-order sum of the pphh ring diagrams 

is H3 

ae™ = ( dx E y ™ (y> A ) y ™ ( fc ^ A ) G £U ^ 

m ijkl(<A) 

where the transition amplitudes Y m and eigenvalues io m 
are given by the following self-consistent RPA equation: 

E[( e i + £ j)%e/ + A(l - m - nj)L ijt ef(u)]Y m (ef, A) 

ef 

= (j, m (u,\)Y m (ij,\); (i,j,e,f)<A. (14) 

The index m denotes states dominated by hole-hole com- 
ponents. The vertex function L is obtained from 2- and 
1-body diagrams first order in G M • The above equa- 
tion is solved with the self-consistent condition that the 
energy variable of L is equal to the eigenvalue, namely 

u = Hm(w, A) = w~(A). (15) 

Comparing with the Vi ow ^u ring diagram calcula- 
tion described earlier, the above G-matrix calculation is 
clearly more complicated. Because of the energy depen- 
dence of the interaction G M , the above equations have to 
be solved self-consistently both for the s.p. spectrum and 
for the RPA equations. To attain this self consistency, 
it is necessary to use iteration methods and this pro- 
cedure is often numerically involved. In contrast, ring- 
diagram calculation using the energy-independent inter- 
action Vi aw —k is indeed much simpler. As mentioned ear- 
lier, we shall carry out ring-diagram calculations using 
both methods. 



III. Vio W -h WITH INFINITE SCATTERING 
LENGTH 

To carry out the above ring-diagram calculation, we 
need the low-momentum potential Vi ow ~k- Since we are 
interested at neutron matter at and near the unitary limit 
(infinite scattering length), we should have Vi ow -kS of 
definite scattering lengths, including ±oo, so that the 
dependence of our results on scattering lengths can be 
investigated. In the present work, we have chosen a two- 
step procedure to construct such potentials so that the re- 
sulting potentials are close to realistic neutron potentials. 
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AE 0-(K) 




an effective low-momentum T-matrix by 
T low - k (p', PlP 2 ) = V^kip' \ P ) 



(ii) 




FIG. 1: pphh ring-diagram summation in the calculation of 
the ground state energy shift. 



We first construct bare potentials V a based on a realis- 
tic nucleon-nucleon potential; these potentials are tuned 
so that they have definite scattering lengths. Renormal- 
ized low-momentum potentials V l ® w _ k are then obtained 
from V a using a renormalization procedure which pre- 
serves the scattering length. 

We start from the high-precision CD-Bonn (29j 
nucleon-nucleon potential. For this potential, the scat- 
tering length of the 1 5'o channel is already fairly large 
(-18.97 /to), and it is found to depend rather sensitively 
on the interaction parameters. Thus by slightly tuning 
the interaction parameters of the CD-Bonn potential, we 
have obtained a family of 1 Sq neutron potentials of def- 
inite scattering lengths. We shall denote them as V a . 
Our tuning procedure will be discussed in section IV(A). 

Recently there have been a number of studies on 
the low-momentum nucleon-nucleon potential Vi ow -k 
[30L l3ll . HI , [33L [H, Viow-k is obtained from a bare 
nucleon-nucleon potential by integrating out the high- 
momentum components, under the restriction that the 
deuteron binding energy and the low-energy phase-shifts 
are preserved. The Vi ow -k obtained from different re- 
alistic potentials ( CD-Bonn [2{|, Argonne [4(| , Ni- 
jmegen [4l[ and Idaho [i^]) all flow to a unique potential 
when the cut-off momentum is lowered to around 2/to . 
The above Vixyw-h is obtained using a T-matrix equiva- 
lence renormalization procedure [3(| HH, HH, Hj, H|| . 
Since this procedure preserves the half-on-shell T-matrix, 
it of course preserves the scattering length. Thus this 
procedure is suitable for constructing V l ® w _ k , the low- 
momentum interaction with definite scattering length. 
Using this procedure, we start from the T-matrix equa- 
tion 



(16) 

where V a is a modified CD-Bonn potential of scattering 
length a. Notice that in the above the intermediate state 
momentum q is integrated from to 00. We then define 



A 2, Vf ow _ k (p ',q)T low _ k (q,p,p 2 ) 
q dq — 



p 2 — q 2 + i0 + 



(17) 



where the intermediate state momentum is integrated 
from to A, the momentum space cut-off. We require 
the above T-matrices to satisfy the condition 

T(p',p,p 2 ) = Ti ow - k (p',P,p 2 ); (p',p) < A. (18) 

The above equations define the effective low momentum 
interaction ■ . The iteration method of Lee-Suzuki- 
Andreozzi [43,l4J] has been used in calculating V^ w _ k 
from the above T-matrix equivalence equations. From 
now on, we shall denote V l a ow _ k simply as Viow-k- 



IV. RESULTS 

A. Low-momentum interactions and scattering 
lengths 

To study neutron matter at the unitary limit, we first 
need a realistic neutron-neutron interaction that would 
lead to a huge 1 So scattering length a s , and a small ef- 
fective range r e . We obtain such interaction by 'tun- 
ing' the meson mass ra a in the usual CD-Bonn potential. 
The exchange of a lighter meson generates a stronger at- 
traction, therefore making the scattering length a s more 
negative until a bound state is formed. As one 'tunes' 
across the bound state, a s will pass from —00 to +00 , 
eventually become less and less positive. In this work, 
this to ct 'tuning' is taken as a manual adjustment in the 
strength of the neutron-neutron potential. Of great in- 
terest is that this 'tuning' may naturally come from the 
density-dependence of the nucleon-nucleon potential via 
the mechanism of Brown- Rho (BR) scaling (45I I46L l47l |. 
which suggests the in-medium meson masses should de- 
crease. 

At normal nuclear matter density, the meson masses 
of , u> and a are all expected to decrease by about 15% 
|47| compared to their masses in free space. This decrease 
will enhance not only the attraction from a but also the 
repulsion from p and w. As a preliminary study, we shall 
tune only m a in the present work. To compensate for the 
repulsive effect from p and u> (which are not tuned in the 
present work), we shall only tune m a slightly, namely a 
few percent. We shall consider that the above BR scaling 
is compatible with neutron matter of moderate density 
(kp ~ l/m -1 ). In a future publication, we plan to carry 
out further studies, including the tuning of p- and to- 
meson masses. 

Various 'tuned' CD-Bonn potentials are listed in Table 
ITTl From there one can see the sensitivity of the scattering 
length to the change in m a . At m a w 442 M eV, namely a 
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name 


m (MpV) 
Ilia I Ivi c V J 


a s {fin) 




original CD-Bonn 


452 


-18.97 


2.82 


CD-Bonn-10 


460 


-9.827 


3.11 


CD-Bonn-42 


447 


-42.52 


2.66 


CD-Bonn-oo 


442.85 


-12070.00 


2.54 


CD-Bonn+oo 


442.80 


+5121.00 


2.54 


CD-Bonn+21 


434 


+21.01 


2.31 



TABLE II: m CT in the original CD-Bonn potential is tuned 
to give neutron-neutron potentials with different scattering 
lengths. 



2.4% decrease from the original, a s w — 12000/m. Notice 
that the effective ranges for the CD-Bonn potentials are 
larger than the actual ranges of them. For example, r e for 
the original CD-Bonn potential is 2.82/m, considerably 
larger than the range of one-pion exchange. Within the 
range of Fermi momenta from 0.8/m _1 to 1.5 fm that 
we use in our computation below, a s ~ —12000/m is 
obviously enormous compared to any length scale in the 
system, thus we expect the neutron matter to be at the 
unitary limit, i.e., no different from the limiting case a s = 
—oo. For convenience, we name such potential CD-Bonn- 
oc. 

Following the renormalization procedures as already 
described in Section III, we obtain the low-momentum 
potential Vi ow -k's for several CD-Bonn potentials listed 
above. A comparison of the diagonal matrix elements 
in the V/ OM ,_fc's (with a fixed cut-off momentum A) is 
shown in Figure O It is of interest that the strength of 
Vi ow -k only changes weakly with the scattering length. 
For example, it changes by merely about 10% from a s = 
-18.97/m to -12070/m. 



0.2 




k/frrf 1 

FIG. 2: Diagonal matrix elements of Vi ow -k constructed from 
CD-Bonn potentials with different scattering lengths. A = 
2.4/m -1 is used in all cases. 



B. Ground-state energy and 



the universal constant £ 

Here we shall present our major results, namely the 
ground state energies E Q of neutron matter at and 
close to the unitary limit from the summation of low- 
momentum ring diagrams to all orders. Following the 
potential renormalization procedure described in section 
III, we first calculate Vi ow -k for certain chosen values 
for the decimation scale A. Then the all-order sum of 
the pphh ring diagrams are calculated using the above 
Viow-k- As introduced in Section II, the calculation de- 
tails in the summation of pphh ring diagrams can be 
found in Ref . [37l] . How to choose the decimation scale 
A is clearly an important step in our calculation, and in 
the present work we shall use a stable-point, or 'fixed- 
point', criterion in deciding A. Before discussing this 
criterion, let us first present some of our results for the 
ground-state energy per particle (Eq/A). In Fig. 3 we 
present such results for four a s values, calculated with 
As determined by the above criterion. (The details of 
this determination will be described a little later.) As 
shown by the figure, we see that Eq/A does not change 
strongly with a s . The ratios £ = Eq/E^^ are then 
readily obtained, as shown in Fig. 4. It is of interest that 
the ratios for the four a s cases are all weakly dependent 
on kp. To help understand this behavior, we plot in 
Fig. 5 the potential energy per particle PE/A (namely 
AEq P I A of Eq.(6)) versus kp, for the same four a s cases. 
It is rather impressive that they all appear to be straight 
lines. We have fitted the 'lines' in the figure to the equa- 
tion PE/A = (h 2 /m) (/3fe|+7): We have found (/3,j) 
= (-0.1370, 0.0002), (-0.1498, -0.0008), (-0.1649, -0.0035) 
and (-0.1797, -0.0082) respectively for a s = -9.87/m, 
-18.97/m, -12070/m and +21.0/m. The rms devia- 
tion for the above fitting are all very small (all less than 
0.0013), confirming that they are indeed very close to 
straight lines. The above results are of interest, and are 
consistent with those shown in Fig. 4. In fact the ratios 
of Fig. 4 are determined by the 'slopes' of these 'lines'. 

Before further discussing our results, let us now address 
the question of how to determine the decimation scale 
A. There are basically two considerations: The first one 
concerns the experimental NN scattering phase shifts on 
which realistic NN potentials are based. The second is 
about the dependence of our results on A. Realistic NN 
potentials [29|, 0, Hjl H2] are constructed to reproduce 
the experimental NN phase shifts up to Ei a i, w SOOMeV. 
This suggests that A is about 2/m -1 , as beyond this scale 
NN potential models are not experimentally constrained 
and are thus rather uncertain (model dependent) [341 ] . 

We now turn to the dependence of our results on A. As 
described in Section II, Vi ow -k is used in the determina- 
tion of the H.F. single particle spectrum (see EqJ5]), the 
transition amplitudes Y in the RPA equation (see EqJTJ), 
and finally, the ground state energy E (see Eq. [6]). In- 
tuitively, Eq should exhibit a non-trivial A-dependence. 
For various Fermi-momenta, this dependence is studied 
and is found to be remarkably mild. 



-16 



FIG. 3: Ground state energy per particle, Eq/A, of neutron 
matter with various tuned CD-Bonn potentials, computed 
from the summation of low-momentum pphh ring diagrams. 
Only So contribution is included. 



0.6 0.1 




0.8 0.9 



FIG. 4: The ratio E /E^ re ' ' as a function of Fermi momen- 
tum kF for the various CD-Bonn potentials listed in Table HT1 
The data with CD-Boon-oo(a s = — 12070/m) indicates that 
E /E^ ree is a constant of 0.443 ± 0.006 over the range of hp 
as shown. 



As an example, let us present in Fig. 6 our re- 
sults obtained with the potential CD-Bonn-oo. For 
A = (2.0 — 2.6)/m _1 , it is seen that £ varies actually 
by a rather small amount (note that the range of our 
plot is from 0.438 to 0.444). Furthermore the A depen- 
dence of £ shows up as a curve with a minimum. The 
final choice of A is based on the criterion that Eq should 
be stable against changes in A. As shown in the fig- 
ure, an obvious stable-point, or fixed-point, defined by 
dE (A)/dA = 0, is found at about 2.3/m _1 . Thus we 
have used A = 2.3/m~ 1 for CD-Bonn-oo. We found that 
the position of the fixed point is almost the same for the 
different Fermi-momenta in the range (0.8 — 1.5)/m . 
The same procedure is done on the original CD-Bonn, 



1.2 1.4 

kp / frrf 



1.6 



2.2 



FIG. 5: Potential energy per particle, PE/A, of neutron mat- 
ter with various tuned CD-Bonn potentials, computed from 
the summation of low-momentum pphh ring diagrams. Only 
1 So contribution is included. 



and other tuned potentials. The fixed points, also with an 
negligible dependence on kp , are found to be 2.15/m~ 1 , 
2.25/m -1 and 2.4/m -1 respectively for CD-Bonn poten- 
tials of scattering lengths — 9.8 fm, — 18.9/m (the orig- 
inal CD-Bonn), and +21.01/m. The above fixed-point 
A's have been used for the results presented in Figs. 3-5. 




2.2 2.3 2.4 

A / flTT 1 



FIG. 6: Determination of the fixed point where dEo/dA = 
for CD-Bonn-oo. 



Of great significance is the ratio of the ground state 
energy to that of the non-interacting case, namely 
Eo/EQ Tee . At the unitary limit, it is expected to be an 
universal constant, named £. This constant is of great 
importance as it determines the equation of state of all 
low-density cold Fermi gas. At the unitary limit, our data 
on E /E* ree all lie within a narrow window from 0.437 to 
0.448. Such result confirms a universality over Fermion 
density in a wide range (1.73 — 11.40) x 10 _2 /m~ 3 . Most 
importantly, the numerical value of £ is remarkably close 
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to that from Monte Carlo methods, which by far is be- 
lieved to be the best estimate. Astra et. al. obtained 
0.42(1) based on a square well potential and particle den- 
sity uRq — 10~ 6 (where Ro is the potential range). Carl- 
son et. al. obtained 0.44(1) based on a 'cosh potential', 
and particle density n/i~ 3 = 0.020 (where 2//J, is the ef- 
fective range). In our case, nA~ 3 = (1.4 — 9.4) x 10~ 3 
(where A = 2.3/m -1 is the decimation scale in the renor- 
malization). These works, including ours, employ very 
different interactions and various particle densities. Still, 
the value of £ agrees incredibly well. 

In Figure [4] we contrast the data from CD-Bonn-oo 
with that from the original CD-Bonn and other tuned 
potentials. Even though the 1 Sq scattering length in 
the original CD-Bonn is already fairly large (a s = 
— 18.97/m) , still the equation of state, as predicted from 
the ratio Eo/EQ ree , has significant difference from the 
unitary limit. As seen in our data with CD-Bonn-oo po- 
tential, at the unitary limit the ratio Eq/Eq 7 "^ — 0.44 is 
practically independent of the underlying neutron den- 
sity n. 



0.55 




0.45 



0.4 



0.35 



0.3 




0.9 



low-k 



■ a =-18.97fm, Ringw/V, 

s » w | 

-a = -18.97 fm, Ring w/ G-matrix 
a = -12070fm, Ringw/V, 

s ° low-k 

-a = -1 2070 fm, Ring w/ G-matrix 



1 



1.1 



1.2 



1.3 



1.4 



k p / fm" 1 



1.5 



FIG. 7: The ratio E Q /E^ ree for the potentials CD-Bonn-oo 
and CD-Bonn(-19.87) computed with two methods. Ring 
w/ G-mat: pphh ring diagrams summation with Brueckner 
G M -matrix. A = 2.3/m~ 1 is used, computation is fully self- 
consistent. Ring w/ Viow-k- pphh ring diagrams summation 
with Viow-k, fixed point is at A = 2.3. 



C. Comparison with G-matrix results 

As discussed in section II, our ring-diagram calcula- 
tions are based on a model space framework. A model- 
space is defined by momentum {k < A} where A is the 
decimation scale. The space with k > A is integrated out, 
resulting in a model-space effective interaction V e ff ■ We 
have used so far the energy- independent Vi ow ^k for V e //. 
Alternatively, on can also use the energy-dependent G M - 
matrix (of section II) as V e f / . These two approaches are 
formally equivalent. We have carried out calculations to 
check this equivalence. 

We have repeated the ring diagram summation with 
the energy-independent Vi ow -u replaced by the energy- 
dependent model-space Brueckner G M -matrix, and carry 
out a fully self-consistent computation in summing up the 
pphh ring diagrams. The exact procedures in Ref.[37j 
are followed (section II). Ring diagrams within a model 
space up to a cut-off momentum A is summed to all or- 
ders. We found that the ground state energy is rather 
insensitive to the choice of A. See Figure [7] for the 
data of CD-Bonn-oo and CD-Bonn(-18.97), done with 
A = 2.3/m~ 1 , 2.25/m~ 1 respectively. As illustrated, 
the two methods, namely, ring diagram summation with 
Viow-k and that with G M -matrix, are fully consistent. 
This is a remarkable and reassuring result, as the cal- 
culational procedures of them are vastly different. For 
the G M case, the s.p. spectrum, the RPA amplitudes Y 
and energies ui^ are all calculated self-consistently, while 
for the Vi ow -k case no such self-consistent procedures are 
needed. Clearly the Vi ow -k ring-diagram method is more 
desirable. 



D. Schematic effective interaction at unitary limit 

At the unitary limit, the simple equation of state 
Eq = iEl ree in neutron matter suggests a very counter- 
intuitive nature in the underlying system: strongly in- 
teracting fermions essentially can be described by a non- 
interacting picture with an effective mass. This unex- 
pected 'simplicity' can best be captured by a schematic 
interaction. To illustrate this, let us consider neutron 
matter confined in a closed Fermi sea |$o(^f))- In other 
words, we consider neutron matter in a one-dimensional 
model space. We denote the effective interaction for this 
model space as V F s- Then the potential energy per par- 
ticle is 



PE 



= (* {k F )\V FS \*o(kf))/A 



where k is the relative momentum. 

Suppose we take Vfs as a contact effective interaction 



V FS = 



1 



(20) 



(h = m = 1) where <S* is a positive parameter with S << 
\a s \. When S=l and k F replaced by A, V F s is the same 
the effective interaction for the pion-less effective field 
theory [34|, [3f|. Substituting the above into Eq.(19) gives 



1 



(21) 



2 a^kr. 
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At the unitary limit (infinite a s ), the above gives 
£=4/9, independent of hp, which is practically the same 
as the result for £(-12070) of Fig. 4. The above also 
gives £ for finite a s . At the unitary limit, we expect 
Vfs to be unique. For finite a s (away from the unitary 
limit), it is not expected to be unique and the param- 
eter S is expected to depend on the underlying poten- 
tial. As shown in Fig. 4, we have calculated £ using the 
CD-Bonn potentials of finite scattering lengths. These 
results can also be qualitatively described by the above 
equation. For instance, for S = 1.25 and k F — 1.0, the 
above equation gives £= 0.54, 0.50 and 0.39 respectively 
for a s =-9.87/m, -18.97/m and +21.01/m. In short, 
certain main features of our results obtained from ring- 
diagram calculations with the CD-Bonn potentials can 
be qualitatively reproduced by the above simple contact 
effective interaction. 



momentum interactions Viow-k of the same specific scat- 
tering lengths. The ground state energy Eq of neutron 
matter are then calculated by summing up the pphh ring 
diagrams to all orders within the model space {k < A}. A 
fixed-point criterion is used to determine the decimation 
scale A. We have carried out ring-diagram calculations 
using two types of renormalized interactions, the energy- 
independent Viow-k and the energy-dependent G-matrix, 
with results given by them being nearly identical. The 
Vi ow -k ring-diagram method has a simpler formalism and 
is also more suitable for numerical calculation. For the 
CD-Bonn-oo potential, the ratio Eo/EQ Tee is found to be 
very near a universal constant of 0.44 over the neutron 
density range (1.73 - 11.40) x 10~ 2 /m~ 3 . Our result 
agrees well with the recent experimental measurement 
and Monte-Carlo computation on cold Fermi gas at the 
unitary limit. 



V. SUMMARY 

In conclusion, we have carried out a detailed study on 
neutron matter at and close to the unitary limit with a 
low-momentum ring diagram approach. By slightly tun- 
ing the realistic CD-Bonn potential, we have obtained 
1 5o neutron potentials of specific scattering lengths, in 
particular the CD-Bonn-oo one with a s of — 12070/to. 
By integrating out their momentum components be- 
yond a decimation scale A, we obtain renormalized low- 
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